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1. Introduction 

Entanglement entropy]^, has become a powerful diagnostic tool for a variety of phe¬ 
nomena in quantum theory (see for example Refs. Very important progress in our 

understanding of entanglement arose from the proposal of holographic entanglement 
entropy, which has beautifully linked the study of entanglement to gravity and thermody¬ 
namics (for example, the recent work of Ref.|^]). 

To define entanglement we break up the degrees of freedom into two disjoint sets A 
and B. The density matrix of the full system is taken to be p. Define the reduced density 
matrix in the region A as pA = tr^p, obtained by tracing out the degrees of freedom in 
region B. Then the entanglement entropy is defined as: 

5^; =-trpA log PA (1.1) 

A more general quantity called the Renyi entropy is often studied. This is defined as: 

Sn = -^logtrp)^ (1.2) 

1 — n 

where n is a positive integer. Besides being a diagnostic on its own, the Renyi entropy can 
sometimes be used as a trick to compute the entanglement entropy. This happens when Sn 
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can be unambiguously continued to arbitrary real values of n. Then it is easily verified that 
the limit lim^^i Sn = Se- The advantage of focusing on Renyi entropy is that it can often 
be computed using the replica trick|^, |l^, H, |T^. For this, one constructs “twist fields” 
and computes their correlation functions using relatively standard procedures in CFT. The 
twist correlators determine a quantity called the “replica partition function” Zn, in terms 
of which the Renyi entropy is: 

‘S'n = :r—— log (1.3) 


, Zn 


Entanglement and Renyi entropies have been computed for a wide variety of quantum 
field theories but they tend to have the problem that either (i) they are very universal, 
depending on limited information about the QFT (in which case they are not a very useful 
diagostic), or (ii) they are extremely hard to compute. In class (i) one has conformal held 
theories at hnite temperature or in a hnite space, with a single entangling interval, for 
which the entanglement entropy depends only on the central charge. However when the 
temperature and spatial size are both hnite, there is a hnite space with multiple entan¬ 
gling intervals, the same theories do not exhibit such a level of universality and therefore 
this situation is more interesting. The case of multiple intervals at zero temperature was 
investigated in Refs.|T^, 14, 15, |T^. At hnite temperature and hnite size (i.e., on a Eu¬ 


clidean torus), only a few explicit computations exist and these apply to free held theories 
of bosons or fermions. The results for the Renyi entropy are found to depend sensitively 
on the spectrum of operator dimensions and not just the central charge. To date, explicit 
computations have been carried out only for the free Dirac fermion ||T^, [l^ and a pair of free 
scalar helds |19, 20, 21|. There are also some general results about the universal thermal 


correction in generic 2d CET’s[| 


]. Quantum corrections to such entropies have been 


addressed in the holographic context in Refs.[§J, |2^. 


Once we are on the torus (dehned by the inverse temperature /3 and the system size 
L), the issue of modular invariance naturally arises. Early computations including those 
in 0,0, obtained the Renyi entropy for free Dirac fermions in a hxed spin structure 
(which could be any of (—, —), (-I-, —), (—, -|-) but not (-1-, +)). The computation made use 
of local Bose-Fermi equivalence to construct the twist held and compute its correlators, 
then manually restricted the bosonic sum to obtain a hxed fermion spin structure. The 
results are not modular-invariant and therefore do not respect Bose-Fermi duality. This 
brings into question whether Renyi and entanglement entropies are properties of a QFT 
independent of its presentation, or depend on the specihc helds used to dehne the QFT. 

This point is discussed in Ref.[^]^ which makes a strong case that entanglement should 
depend only on the theory rather than on its presentation. The consistent dehnition of 


^We thank Shiraz Minwalla for bringing this work to our attention. 
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a CFT incorporates modular invariance as a fundamental property and this provides a 
valuable constraint on its operator content, partition function and correlators [26, 


Motivated by this viewpoint we revisit the previous computation of Renyi entropy for 
free fermions, performed using the replica trick, and examine whether a modular-invariant 
replica partition function exists satisfying all the desired properties. In order to achieve 
modular invariance one needs to perform a sum over spin structures (i.e. over anti-periodic 
and periodic boundary conditions in both the spatial and thermal directions) 


The replica trick involves extending the original theory to an n-fold copy and per¬ 
forming a twist-field computation on this copy. The presence of fermion spin structures 
introduces, quite literally, a new twist on the problem - one has to decide whether the spin 
structures are correlated or uncorrelated across replicas. To get some physical input into 
which one is correct, we appeal to the universal relation between entanglement and thermal 
entropy first noted in Ref.jj^. We find a perhaps surprising result: for free fermions, this 
relation holds only if we choose two different ways of summing in two different limits: for 
a small entangling region one must consider uncorrelated spin structures, while for a large 
one (close to the size of the entire space) one has to completely correlate spin structures 
across replicas. The correct formula for a general size of the entangling region must in¬ 
terpolate between these two extremes, and we speculate at the end about how this could 
come about. 


It is known that for large intervals, the replica trick requires a correlation between 
states on different replicas, as was originally evident from the work of p^ . In a beautiful 
series of works, this relation has been developed further and a formal proof been provided 


for the entanglement-thermal entropy relation 0 


]. Our computation can be seen as a 


verification of these ideas in the specific context of fermion spin structures. In Ref.[21] the 
proof has been explicitly verified for the free-boson partition function. We confirm that our 
result agrees with that of the above reference in the special case of radius R = 1 (in units 
where a' = 2), where the free boson is known to be equivalent to a free Dirac fermion. 
This provides strong evidence, for the first time, that Bose-Fermi duality is respected by 
entanglement entropy on the torus. 


Both the ways of summing over spin structures that we use lead to a replica partition 
function that is modular covariant rather than invariant. Under modular transformations of 
the torus, it transforms into itself with a prefactor. As a result the Renyi and entanglement 
entropies acquire an additive contribution upon modular transformations. We will see that 
the covariance arises from the fact that twist field correlators are involved: correlation 
functions in CFT are modular covariant rather than invariant, and the precise prefactor 
that we find agrees with the general expectation. We verify that the Renyi entropy for 
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a compact free boson at radius R, originally attempted in []^| and recently corrected in 
Ref.( 23 ], is also modular covariant with the expected prefactor. 

Next we turn our attention to more general free 2d CFT’s. The CFT of two Dirac 
fermions with correlated spin structures is equivalent to a pair of bosons compactified at ra¬ 
dius R = \/2 (this is the self-dual value under T-duality). We compute the Renyi entropies 
for this theory for small and large intervals and find modular-covariant answers that satisfy 
the thermal entropy relation. This is then repeated for theories of multiple correlated Dirac 
fermions. These are CFT’s of arbitrarily large central charge d that are not direct sums of 
simpler CFT’s. The bosonic duals are given by multiple free bosons compactified on the 
weight lattice of Spin(2d) with a specific background R-field|^. We use the weight-lattice 
structure to construct a twist field, and compute its two-point function. Incorporating the 
sum over spin structures in the two relevant limits provides the replica partition function in 
these limits and confirms the arguments made above. We conclude with some observations 
and speculations about future directions. 


2. Free Dirac fermion and modular invariance 

Consider a CFT whose partition function on a rectangular torus of spatial size L and 
Euclidean time extent /3 is Zi(L,/3). As mentioned above, the Renyi entropy Sn associated 
to a spatial region A running from 0 to £ is defined as: 

S'n = log tr pd (2.1) 

where pA is the density matrix obtained by tracing out the degrees of freedom outside 
A. This can be evaluated using the replica trick, which we briefly summarise here. One 
extends the original torus to an n-fold cover with branch cuts along spatial intervals from 
0 to i. This can equivalently be thought of as a single copy of the original torus with the 
insertion of “twist fields” labelled at the endpoints of the entangling interval. Here k 
ranges from to in integral steps. The job of the twist fields is to introduce a 

phase factor e~^ on the free fermion field, which we denote D{z), as it goes around the 
twist field: 

ak{z, z)D{w) ^ {z — w)^ (2.2) 

and the opposite phase on the anti-holomorphic conjugate fermion: 

ak{z,z)D{w) ^ (z - w)~« (2.3) 

The twist fields can be shown to satisfy = -^ {it- — ^) hi any CFT, where c is the 

central charge. 
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One can then show that[^, |l 2 t 


txp\ 


n — 1 
2 


n 


C7_fc(0,0)) 


k=— 


n — 1 
2 


(2.4) 


It is convenient to think of the product of un-normalised correlators as defining the “replica 
partition function” 


Zn. — 


n Zi{ak{i,£)(T_ki0,0)} = (((Tfc(£,£)(T_fc(0,0))) 


(2.5) 


fc =-2 


k=-^ 


where Z\ is the ordinary partition function of the theory (without insertions), and we use 
((• • •)) to denote un-normalised correlators. It follows that: 

= ^ ( 2 - 6 ) 

from which the Renyi entropies are easily obtained. 

Consider a Dirac fermion with a single entangling interval of length i. This theory con¬ 
sists of two Majorana fermions with correlated spin structures. Labelling the holomorphic 
parts of the Majorana fermions as V’l) ^^2 we make the Dirac fermion D{z) = 'il>i{z) + i^p 2 iz). 
The local operators of dimension (^, 5 ) arising from this are: D{z)D{z), D^{z)D{z), 
D{z)D^(z), D^{z)D^{z). This theory has the modular-invariant partition function: 


Zbi: 


E 

i/=2,3,4 


r]{T) 


(2.7) 


For future use we define: 


Zi[m] 



u=2,3,4: 


0u{0\t) 

r]{T) 


( 2 . 8 ) 


where the subscript 1 refers to the ordinary partition function (this will later be replaced 
by n for the reth replica partition function). This is the ordinary partition function for m 
Dirac fermions with correlated spin structures. In this notation the single Dirac fermion 
has ordinary partition function Zi[l]. 

For the replica partition function of the same theory, the following result^ was obtained 
in Ref, : 


((Tk(z,z)a-k(0,0)} 


0j(O|r) 

4Afe 


0 i(f|r) 


9u{0\t) 


(2.9) 


^In their conventions, L = 1 and i is denoted by L. 
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Here A^. = i/ = 2, 3,4 is a fixed spin-structure on the torus corresponding to boundary 
conditions (-I-, —), (—, —), (—, -|-) respectively^, and t = For the (-I-, -|-) spin-structure 
the result would formally be the same but with u = 1. Subsequently the above results 
were re-derived in Ref, , where it was observed that the 2-point function in the (-I-, -|-) 
spin structure is divergent due to the vanishing of 0i(O|r). Importantly, the entanglement 
entropy following from the computations in Refs. MM satisfies the relation proposed in 
Ref. 117], based on holography, relating its small and large-interval limits to the thermal 
entropy of the same system. 

Despite meeting all the other physical requirements, the above result raises a puzzle. 
As emphasised in Ref.[^], entanglement entropy should be a feature of a definite quantum 
field theory and independent of the presentation of that theory. If so, in a theory where 
Bose-Fermi correspondence holds, one should obtain the same result in both bosonic and 
fermionic formulations. A free fermion with a single spin structure is not dual to the 
theory of free bosons, rather it is the modular-invariant free fermion partition function 
(and correlation functions) that can be compared with those of the boson theory. In this 
spirit, we propose that Renyi entropy be computed by performing a sum over all four spin 
structures in the replica partition function and dividing by a corresponding quantity 
in the absence of replicas. In principle there is more than one way to do this, and this 
point will be central to our discussion. 


Let us review the calculation of [lAj for the Dirac fermon theory. To compute the 
correlators of fermionic twist fields on the torus in this theory, these operators are first 
identified with (non-local) operators in a free boson theory using the fact that a Dirac 
fermion is equivalent to a free compact boson of radius R = 1 (in our conventions a' = 2, 
so T-duality acts by R —and R = \/2 is the self-dual radius). At a general radius R, 
the free compact boson has vertex operators labelled by integers (e,m): 


de,m(^)2) — Fe,m('2)Ve,m(^) 


( 2 . 10 ) 


where 


of conformal dimension: 


V„(^) = e'(*+")*W 

V„(z) = 


( 2 . 11 ) 


— \2\R^ 2 ) ^ 2\R 2 ) 


( 2 . 12 ) 


®The first entry labels the spatial boundary condition and the second, the boundary condition in Eu¬ 
clidean time. 
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The OPE’s between the chiral parts of the vertex operators are given by: 


( \ ( rn_B\ 

^ ^ ^Ve+e' 


ra+m 


iz) 


(2.13) 


For the special value i? = 1, the compact boson is equivalent to a Dirac fermion. The 


fermion held has (A, A) = (^ , ^ ) and is given, in bosonic language, by Oi,o- In Ref.[]T| 


the fermionic twist helds are identihed as: 


Ufc = Og 2. (2.14) 

’ 71 

for k = —These operators have (A, A) = (^j^)- Since the winding 
number m for these operators is not an integer, they are not included in the set of local 
operators of the theory. This is expected, since twist operators are mutually non-local with 
the physical helds of the theory. One has the OPE’s: 

k — k 

Og ^{z, z)Vifi{w) ^ {z — w)ri, Oq m{z, z)Vifi{w) ^ {z — w) ri (2-15) 

as desired. Now we only need to compute: 


((Oo^2fc(z,z)Oo,_2io(0))) 


0 ,- 


(2.16) 


For this we use the general result (see for example Ref.||^): 

4 Ae. 7 Ti 


{{Oe,m{z.,z)0-e-m{^))) = 


e'Mr) 


1 


\'n{r)? frt, 


g2A^, -2A^, g47ri(a 


e', OL^i 


(2.17) 


With (e,m) = (0, —) and ii = 1, this reduces to: 


i0^i_k{z,z)0^_2_k{^))) = 


O'Mr) 

71^ 1 

0M\r) 

r/(r) 2 ^ 

' ' ^ ^' e,m 


2fe2 

9Mr) 

^ ^ 1 ELi 


^i(lk) 




kl 


(2.18) 


An important point to emphasise here is that in this approach one has to use free bosons 
to represent fermionic twist helds. These do not introduce a phase on the boson, instead 
they shift the free boson (f appearing intp = by sending ^ which is the /c/n’th 

multiple of its natural period 27r (because iZ = 1). Correspondingly one has 4> ^ 

(here are the holomorphic and anti-holomorphic parts of the scalar held respectively). 
The bosonic representation is used only as a tool to obtain the correlators of the fermion 
twist helds. 
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At this point, we need to decide how to take the product over replicas. The most 
straightforward way would be to just take the product of the above result over all k. In 
this way the spin structures are summed over before we carry out replication. Hence the 
spin structure on one replica of the torus is uncorrelated with that on any other, leading 
to what we call the “uncorrelated replica partition function” Z^: 


n—l 

2 


’ n ’ n 

l(n-i) 

6 ^ n' 2 


k=-'- 


(2.19) 


0 'i(O|r) 




n 


k=-'- 


2 |,7(r)|2 


There is another way that one can take the product. First isolate the contribution to 
Eq. ( p.l8 ) corresponding to a given spin structure v G 1,2,3,4. This simply corresponds 
to picking out the function 9^. Next, take the product of this term over all replicas and 
finally sum over spin structures. In this way the spin structure on each replica of the torus 
is the same. We call this the “correlated replica partition function” 


Z^(L,/3;£) = 


1 


0j(O|r) 




1 1 \ n — l 

6 n 1 4 2 

E n 




!2=1 1 - ri -1 

2 


|r/(r)p 


( 2 . 20 ) 


( 2 . 21 ) 


Notice that the two types of replica partition functions coincide at n = 1: 

..c iEh,|9.(0k)P 

Mryr~ 

which is the ordinary modular-invariant partition function of a Dirac fermion (the contri¬ 
bution from v = 1 actually vanishes in this case). We also observe that as 1" ^ 0 (and 
generic n), the two types of partition function are quite different: 




Z;;(L,/3;£^0) 


[2 |r?(r)|2 

1 Y:t=iK{o\Tr^ 


( 2 . 22 ) 


2 |^(e)P" 

The second factors in the two cases are the ordinary partition functions of n Dirac fermions 


with, respectively, uncorrelated and correlated spin structures!291. 

In principle we can consider several more variations in which we correlate and sum over 
the spin structures on various subsets of the n replicas, and let them remain uncorrelated 
across replicas. Such quantities are also modular covariant as one can verify. For the 
moment we restrict ourselves to these two “extreme” cases and will comment about the 
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intermediate cases later on. We could also consider variations in which we sum over a subset 
of the spin structures in a correlated fashion and the remaining ones in an uncorrelated 
fashion. However, for modular covariance it is essential to sum over all spin structures in a 
symmetric way. Hence this possibility is not covariant and we will not consider it further. 

Returning to the two replica partitions described above, we can define two Renyi 
entropies: 


5“ = lim 


1 


5" = lim 


1 — n 

1 


log 






1 — n 


yc 

log- " 




(2.23) 


The denominators are the same because, as we pointed out earlier, the two types of partition 
functions coincide at n = 1. 

Clearly we would like to know which one of these (if any) is the correct Renyi entropy 
of the modular-invariant free Dirac fermion theory. Before addressing this question, which 
turns out to be quite subtle, let us establish a few properties of the two quantities and 


Z'^. First of all, they can be compared to the non-modular-invariant result of Ref.|^| for 
the replica partition function; 


= 


m\r) 


0i(f|r) 


Hn-i) ^ 

6 ^ n' 2 

n 




h— n — 1 


|r?(r)p 


(2.24) 


and the corresponding Renyi entropy: 

5(C = 


1 — n 


log 


7(^1 

(C'.i 


(Zl 


(2.25) 


( 1 ) 


This makes sense only for the spin structures 2,3,4, since Z\ vanishes and the Renyi 
entropies in this spin structure are correspondingly ill-defined. Notice that the expressions 


in Eq. (2.23) are not made from linear combinations of the ratios in Eq. (2.25). 

Let us now check the modular transformation properties of the expressions in Eq. ( 2.1(]| ) 
and Eq. ( 2.20| ). This transformation exchanges the two cycles of the torus. We also have 
a branch cut along the horizontal axis, which under this transformation becomes a branch 
cut along the vertical axis. Thus the transformation acts as /3 -H- L and (. —>■ i^, where we 
have used the identification r = ir 2 = ij- and z = j. This permits us to use: 




= (-i)“^(-ir)5e^ 


7j3a 


Z,T 


(2.26) 
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with the usual dictionary: On —)■ —0i, ^*10 ^ ^2, ^00 ^ 3)^01 ^ 4 - Applying this to the 

second term in Eq. ( 2.19| ) or Eq. ( 2.20[ ), one finds that even though the two expressions are 
different, they pick up the same multiplicative factor; 


g6r ' 


(2.27) 


Erom the first term, which is common to both expressions, we get a corresponding factor: 


which cancels the previous one, as well as a factor: 

X 6 \ n / 

Thus at the end, one has: 






(2.28) 


(2.29) 


(2.30) 


We see that even after summing over spin structures, the replica partition functions are not 
modular invariant, but rather modular covariant - they acquire a multiplicative pre-factor. 
This factor vanishes at n = 1, so Zi is indeed modular invariant as it must be. 

The origin of the pre-factor lies in the fact that the replica partition function Zn is not 
a partition function at all, but rather a correlation function - of twist fields. On the torus, 
the two-point function of a field (j){z, z) of conformal dimension (A, A) transforms as (for 
this and other standard results in GET, see for example Ref. |pCI|| ): 

(2.31) 


cr-\-d ’ CT-\-d 




ar + b 

+ d 


Erom this it is easy to verify that the product of twist-field correlators should transform 
precisely as we have found above. Indeed, since the dimensions (Afc,Afc) of twist fields 
satisfy the universal relation: 

_ n f 1 \ 

(2.32) 


1 

n- 

n 


VAfc = — 

24 

k 

one expects that for every GET of central charge c the analogous result will hold with: 

'j3\ 




L 


Zn{L,p-e) 


(2.33) 


In what follows, we will explicitly verify this in all known cases. 

One still has to physically understand the impact of this modular covariance on the 
Renyi entropy Sn- Under Eq. (^.33 ), one finds that Sn that transforms as: 


S4f3, L) = log ^ + Sn{L, /3) 

bn L 


(2.34) 
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The additive term survives in the limit n ^ 1, therefore it even appears in the entanglement 
entropy. Now we know that the entanglement entropy always contains a non-universal 


additive constant (see for example Ref.|]^]), and one can take the point of view that 
this constant simply changes under modular transformations. Such an interpretation is 
supported by the fact that the additive term is independent of the length I of the entangling 
interval. 

Alternatively, we can modify the prescription for computing the replica partition func¬ 
tions (in both uncorrelated and correlated cases) by multiplying them by an overall factor 
that renders them modular invariant. This is done as follows (importantly the factor is 
independent of £): 


z^ = \j 


Zn 


(2.35) 


We have put in a dependence on c, the central charge, which is equal to 1 in the present 
example. One expects that the replica partition function for every CFT, after multipli¬ 
cation by the above factor, will be modular invariant. From the discussion above, such a 
modification of Zn can be achieved by re-normalising all the twist fields as: 




O'k 


(2.36) 


using the fact that — -). 

Returning now to the Dirac fermion, the two Renyi entropies, defined as: 


QU,C _ 


1 


log- 


(2.37) 


(2.38) 

1 in the above will also 


1 - n {ZiY 

will satisfy: 

S))-^(/3,L;z^) = S"’^(L,^;^) 

It follows that the entanglement entropies obtained by taking n 
be modular invariant. 

In the following subsection we address the question of which one of S')), S')) is the correct 
Renyi entropy of the theory. From now on, we will not be very careful to insert the prefactor 
above or to distinguish between modular covariant and modular invariant expressions at 
each stage, since this insertion is ^-independent and can always be carried out at the end. 


3. Relation to thermodynamic entropy 


In Ref.|17|, the result Eq. (2.24) was used to provide evidence for the conjecture, based on 
holography, that the entanglement entropy satisfies: 


Im {Se{L -i)- Se{Y) = Sth 


(3.1) 
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where S is the thermodynamic entropy of the theory: 


d 




1 


(3.2) 


We will use this as a guide to determine which of our modular-invariant definitions of 
the replica partition function, or Z^, is correct. As emphasised in Refs.p^, |^], we 
can simply use the low-temperature expansion for this purpose. Note that the overall 
factor that we removed from the replica partition function to render it modular-invariant 


is independent of the interval size i and therefore drops out of Eq. (3.1). 

Consider now the small-1" behaviour of the candidate Renyi entropies. From Eq. ( 2.22| ) 
we see that in this limit, —)> {Zi)^ but Z^ (-^i)*^ (upto an overall power of €). It is a 

physical requirement that for small intervals (equivalent to a very large spatial size L) the 
replica partition function should indeed tend to (Zi)^ times the given power of i. On this 
criterion, the correct Renyi entropy for small £ should be 5“ and not S^. 

Next consider the large interval limit, i ^ L. In this case, it has been predicted on 
general grounds |^, that the replica partition function should tend to Zi(nr), i.e. the 
ordinary partition function on a torus with n times the original modular parameter, apart 
from the same power of i as in the ^ ^ 0 case. Moreover it has been shown that this 
behaviour (together with the correct £ —)• 0 behaviour) ensures that Eq. ( |3.1D holds. Let 
us subject our Z^ and Z^ to this test. We start with the latter. From Eq. (|2.20|) one sees 
that: 


1 


= 2 I 




n — 1 
2 

n 


|77(r)|2 

We may now apply an identity involving 0-functions after fractional shifts 


E 

U=1 J - n-1 

2 




(3.3) 


n—1 

2 

n 


9y{- - Z 

n 


n 

.p=i 


(1 - q^py 


1 _ q2pn 


\0i,{nz\nT)\ 


(3.4) 


This is easily derived using the product representation of 0-functions. It is important that 
the above relation is true for the three spin structures i/ = 2, 3,4 with no relative factors. 
The fourth spin structure n = 1 does not matter because we have to put z = 0 in the above 
equation, and 0i(O|r) = 0. Using the above relation, we find: 




1 |0i2(O|nr)p 

2[lJ ^ |r?(nr)P 





Zi(nr) 


(3.5) 
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Thus we see that for a large interval, becomes equal to Z\ {nr) upto the standard power 
of i, as desired. However, if we repeat the analysis for Z^ then this is not the case. With 
this function we would find cross-terms between 0-functions of different spin structures, for 
which an identity like Eq. (|3.4D does not hold. We thus have clear evidence that as i ^ L, 
the correct replica partition function should be Z^ and not Z^. 

Hence we are led to propose that the correct modular-invariant partition function for 
the free Dirac fermion theory is Z^ at small i and Z^ at large i. Let us now verify that 
this proposal precisely satisfies Eq. (p.ip. We have: 


Im (^Se{L -£)- SEii)) 


lim-log 

n^l 1 — n 



Sth 


/ ^i(nT) \ 

\iZi{r)r) 

^z{i) 


(3.6) 


Thus our proposal satisfies the relation to thermodynamic entropy. Notably, the final 
steps work out just as in the general proof of Ref.|2^ (see their Eq.(ll)) but here we 
have used explicit 0-function identities to obtain them. The fact that spin structures are 
uncorrelated/correlated for small/large intervals is a nice mathematical verification, in a 
special class of models, of the general principles put forward in Refs, 
conclusion for us is that neither 5“ nor can be the correct Renyi entropy for all values 
of This naturally raises the question of what is the general answer that interpolates 
between these two limits, a point to which we will return at the end. 

As noted above, in addition to the fully correlated and fully uncorrelated replica parti¬ 
tion functions, one can write down other candidate quantities in which the spin structures 
are correlated within subsets of the total set of replicas and uncorrelated among different 
subsets. It is straightforward to check that such quantities do not behave correctly in 
both the limits considered here, in other words as .£ —>• 0 they do not give {Zi)'^ and as 
i ^ L they do not give Zi{nl3). They are therefore not useful candidates to describe the 
limiting behaviour of the Renyi/entanglement entropies. However, as we will suggest in 
the concluding section, they could play a role for intermediate values of i. 
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The key 


4. Compact bosons 

The replica trick can be applied for free bosons as well. Recall that the partition function 
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for a single free boson at radius R is: 


/ . 1 V—> / e I / e mR 


|r/(r)|‘ 


e,m 


E -^e ^^ m H 

q-W'^q 2 
e m 


(4.1) 


where in the second line we have specialised to real q = . 

For the replica partition function, one considers twist fields Tk with k = 0,1, ■ 
satisfying: 

k — k 

Tkiz,z)(j){w) ~ (z - w)^, Tk{z,z)<p{iu) ~ (z - w)~^ 


and one has: 


n—l 


Zn = Y\{{Tkiz, z)T-k{0,0))) 


, n — 1 

(4.2) 

(4.3) 


k=0 


An important part of the computation of this quantity was carried out in Ref. [p)9|] for a 
pair of free bosons compactified on a square torus of size R. The result is the product 
of a quantum and a classical part. It is rather implicit, involving integrals of products of 
fractional powers of 0-functions which appear in the construction of cut differentials on the 
torus. Unfortunately the classical part of their answer, which carries all the R-dependence, 
is not invariant under R —^ ^ so it does not satisfy T-duality. In Ref.|21] this classical 
part has been corrected and takes a different form that possesses manifest T-duality. 

Accordingly we will analyse the replica partition function obtained in the latter refer¬ 
ence and compare it with our proposal regarding modular invariance. We adapt it to our 
notation, which includes setting a' = 2, and to the case of a single free boson. The latter 
just involves taking the square root of their answer, which is particularly simple since it is 
presented as a perfect square. We write the answer as a product of four factors: 


Z„(R) = ZWZ^)Z^)(R)Z® 


R 


(4.4) 


Only the last two factors, which form the classical part of the answer, depend on the 
compactification radius R. T-duality invariance is already manifest at this stage. Here, 


Z^^'> = 


= 


1 


n—l 

n 


1 


|0(r)|2- 11 |VFii(A:,n;f|r)| 


0j(O|r) 


k=0 


Oliilr) 

Z^\r) = E exp 


(4.5) 


vrR' 


n—l 


2n ^ 
k=0 


Wi{k,n) 

n—l 

Wlik,n) 

j,j'=0 ^ 


cos27r(j -/)- 
n 


rrijiTiji 
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Here and W 2 are integrals of the cut differentials over the different periods of the torus: 


Wlik,n;j;\T)= [\z9iiz\T)-(^-^)ei{z-j;\T) hi {z - ^\r) 

^ (4.6) 

Wiik, n- f |r) = r dz h{-z\r)-hi {z - f h [z - {l - |) f |r) 

JO 

The answer Eq. (ff.4D is normalised such that Zi is indeed the usual partition function of 
the theory. To check this, observe that for n = 1,A: = 0 one has Wi = IjTTf ~ —'iT 2 . 
Hence: 


zW 


1 

ww 


= 1 


z^^\r) = E exp 


ttR^ 


-T2m = 




^(3) 



J^exp 




(4.7) 


Agreement with Eq. ( [4.1] ) is immediate. 

We would now like to investigate the modular transformation of Eq. (4.4). To this end, 
we note the following results: 


r/( - i) = (-iT)2 77(r) 

I - 7 ) =h~^^^^~^'>W^{k,n-,{\T) (4,8) 

_ • -i21^G'l{^\'T) 


The first and third of these are well-known, while the second was shown in Eq.(B.41) of 
Ref. (unlike the very subtle issue of the g-expansion of Wf, this result is straightforward 
to derive). Erom the above, it follows that under a modular transformation, 


W^{k,n) 


Wl{k,n) 

Wl{k,n) 


Wi{k,n) 


(4.9) 


Eollowing this with a multi-variable Poisson resummation as in Ref. j21], we hnd that: 





zi^) 




n 

22 

Rd^ 

n 

22 




Wi{k,n) 2 
Wl{k,n) 

W2{k,n) 2 

Wl{k,n) 



Z® (^R- z 


(4.10) 
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Thus the product transforms as: 


z^\r) 


/n-l 

n 

\k=0 


Wiik,n) 


Wlik,n) 




(4.11) 


(4.12) 


Putting everything together, we find that: 

=|r|H-^)z„(i?;z|r) 

Thus, it is modular covariant with the expected prefactor (see Eq. 

The general expression for Renyi entropy Eq. (4.4) is rather complicated and it has 
proved difficult to extract explicit results from it. However, in Ref.|^^ the limits as ^ 0 
and £ ^ L have been extracted and the results correspond (upto the usual powers of i) to 
(Zi(r))"' and Zi(nr) respectively, where by Zi we mean the expression in Eq. ( [4.11) . This 
of course depends on R, and for R = 1 it coincides with the free Dirac fermion partition 
function. R follows from this and our previous analysis that our fermion replica partition 
function, in the regimes where we know it (small and large i) coincides precisely with that 
of the boson theory at R = 1. R is amusing that the boson theory “knows” about the 
switchover from uncorrelated to correlated replicas within the single, though complicated. 


function in Eq. (4.4). We expect there to be a similar fermion partition function that 
interpolates between the small and large-£ regimes. It should somehow involve an admixture 
of correlated/uncorrelated replicas, as we discuss in the last section. 


5. Other free-fermion CFT’s 

It is interesting to try and extend the previous calculation to free-field CET’s that are not 
direct sums of the Dirac GET. Recall that in the above calculation we had to use both the 
fermionic formulation of the Dirac GET (to identify the right free field to replicate) and also 
the bosonic formulation (to identify the twist field and compute its correlation function). 
Eor multiple Dirac fermions with correlated spin structures, the partition function is not 
simply a power of the single-Dirac partition function^. Gorrespondingly, the bosonic dual 
is not made up of copies of the R = 1 free boson. Indeed in bosonic language it is not the 
product of any number of independent bosons but rather a set of bosons compactified on 
the weight lattice of the algebra Spin(2d) with a specific background R-field turned on|p^. 
We start with the relatively simple case of Spin(4)~ SU(2)x SU(2) and then go on to a 

“^In this part of the discussion there is some potential for confusion: the multiple fermions can have 
correlated/uncorrelated spin structures, which is part of the definition of the theory, while the replica 
partition function can also have correlated/uncorrelated spin structures across replicas as was discussed in 
previous sections. 
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general Spin(2d) lattice. Calculation of the R&yi entropy will again require use of both 
the fermion and boson descriptions of the theory. In each case, the twist field must be 
constructed using lattice vertex operators and their correlation functions computed in the 
theory of lattice bosons. We will again find answers for the replica partition function at 
small and large £ and verify their modular properties. 


5.1 Two correlated Dirac fermions 


= ^ E 


(5.1) 


The theory of two Dirac fermions with correlated spin structures has the ordinary partition 
function!^: 

|0..(0|t)|-‘ 

u=2,3,4: ' 

We would like to find its replica partition function. Let us examine the fundamental fields 
and physical operators in this theory. We have two Dirac fermions that, following the 
notation of Ref.|^], we call Di{z) and D 2 {z). By definition they possess hermitian conju¬ 
gates d\{z), i = 1,2, that are also holomorphic. Combining with their anti-holomorphic 
counterparts, we have 16 physical operators of dimension (^ , ^ ), namely: 


DiD.,D]D.,D,D\,D]Dl 

^ J ’ I J ^ ^ .7 


(5.2) 


where i.j run independently over 1,2 (had the spin structures of Di, D 2 been uncorrelated 
we would only have had the subset with i = j). This c = 2 theory can be bosonised into two 
compact bosons (j)i , (j )2 that are orthogonal to each other and compactified at the self-dual 
radius, which in our conventions is i? = y/2. Thus their periodicity is: 


(pi —y (pi ‘I'sPIk 


(5.3) 


If we denote the two left-moving bosons as pi{z) and the right-moving ones as pi{z), then 
the allowed vertex operators in this theory are: 




(5.4) 


where the indices j,k are independently summed over 1,2. In terms of these bosons, the 
free Dirac helds are represented as: 

Di{z) = D2 {z) = (5.5) 


The Hermitian conjugates have the same expressions but with i —?■ —i, while the anti- 
holomorphic conjugates have pi pi without any change of sign on i. The 16 local oper¬ 
ators obtained after combining left- and right-movers correspond to all pairs of orthogonal 
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vectors of (length)^ = 2 in the unit 2d square lattice: 


(ei;mO = ±((1,1);(0,0)), ± ((1,-1); (0,0)), ± ((0,0); (1,1)), ± ((0,0); (1,-1)), 

±((1,0); (0,1)), ±((1,0); (0,-1)), ± ((0,1); (1,0)), ± ((0,1); (-1,0)) 

(5.6) 

In order to compute the replica partition function we need to find the twist fields which 
satisfy Eq. (|2.2|) . On inspecting Eq. (^), we see that the monodromy induced by the twist 
field is: 


-^Sl -^<*2 


-^4>2 

e ™ e 




(5.7) 


Thus it must leave the boson cj )2 inert, while it shifts (j)i by a fraction ^ of its period: 

2y/2'Kk 


(t)i^ <t)i + 


n 


(5.8) 


We see that the twist field is: 


(Tfc = e e " 


(5.9) 


In other words, it corresponds to . Remarkably this is the same operator that plays 

’ n 

the role of the twist field in the c = 1 theory of a single Dirac fermion - however, in the 
present case it is evaluated at i? = \/2 rather than R = 1. Due to the change of radius, 
it has conformal dimension Ag 2 k = ^, which is twice that in the single-Dirac case. This 
ensures that Ylk ^ “ n) c = 2, as required. 

To calculate its correlator we may use Eq. ( ^.17 ) and make the appropriate change in 
radius, to get: 


{{o.2Az,z)o._2Am = 


e'Mr) 


1 


|r/(r)|^ 


'y ^ g2(Aej^ ,mj+'^e2,m2 ) g2(Aej^ -|-Ae2,m2) 


ei,mi;e2,m2 


e'Mr) 


ei{i\r) 




\r]{T)Y 


The last line follows by the change of variables: 


^ (^^1 ± n2 ± ns ± 77.4), ”71 = ^ (ni ± n2 - ns - 774) 

62 = ^ (ni - n2 ± ns - n4), "^2 = ^ (ni - 772 - ns ± 774) 


(5.10) 


(5.11) 
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The rii can be all integers or all half-integers. Moreover, the sum of all the n* is constrained 
to be even. Upon implementing these two facts, the sum in Eq. ( 5.1C| ) splits into four parts, 
one corresponding to each spin structure. In this way we recover the last line of Eq. ( 5.1C ). 

We now have to take the product over replicas. As before, we can do this after or 
before summing over spin structures. The results are then: 






o'lm 

Oliilr) 



2 |ry(r)|4 


Oiiilr) 



4 


n — 1 
2 


E n 


U=1 1 - n-1 

2 


|ry(r)|4 


(5.12) 


We again propose that the first result gives the correct Renyi entropy as 0 while the 
second one is correct when i ^ L. 

As emphasised at the beginning of this section, the above partition functions do not 
directly follow from the single Dirac case. The reason is that one needs a bosonic dual theory 
in order to write the twist held. For two Dirac fermions with correlated spin structures, 
the bosonic dual is two copies of the R = \/2 compact boson. Interestingly the twist held 
that we found was not symmetric in the two bosons, depending only on one of them. We 
will see below that this pattern repeats for multiple correlated fermions, dual to bosons on 
Spin(2d) weight lattices. 

It is straightforward to demonstrate that under modular transformations the above 
replica partition functions satisfy: 


zr{l3,L) = 




L 




(5.13) 


which precisely agrees with our expectation. 


5.2 Multiple correlated fermions and lattice bosons 

For d > 3, the theory of d free Dirac fermions with correlated spin structures is dual to a 
specihc compactihcation of d free bosons on a target-space torus: 

RC ^ 

where T^ is the root lattice of Spin(2fi) (this can be achieved by starting with a rectangular 
torus and choosing a suitable constant metric and R-held). In this case the d different 
bosons are not orthogonal to each other, while the fermions have correlated spin structures, 
so on both sides of the Bose-Fermi duality we are dealing with CFT’s that are not direct 
sums of simpler ones. 
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The ordinary partition function is well-known and has been derived directly within the 
free fermion representation as well as from the free boson picture|29|. It corresponds to 
Z\ [d] in terms of the definition in Eq. . We would like to construct the replica partition 
functions for these theories. As in the previous examples, this will be possible only by using 
features of both the fermionic and bosonic descriptions together. 

The first step in doing this is to identify the ^ ) operators in terms of the bosons. 
As in the case of two Dirac fermions, we have the operators: 


DpDg, D\Dq, DpD\, dId\ 


(5.15) 


where p,q = 1,2, - ■■ d. In the free boson theory, let A^j be the root lattice and Aw be the 
dual weight lattice. Then the vertex operators are: 


O 






(5.16) 


where w^,w^ G Aw and & Ar. Elements of the weight lattice can be parametrised 

as: 

(5.17) 

where Vi,Vi are integers and is the inverse of gij which is the half the Cartan matrix of 
Spin(2d). Eor the difference to lie in the root lattice, we must require that — Vi) = 

\/2ni where Ui are integers. 

The 2-point function of scalar fields is: 






z)(j){z', z') = -gij log jz - z'f (5.18) 

Hence the conformal dimension of the above operators is (A^,i, A^i) = ^ {gijW^wA gijw’‘w^). 
Thus to reconstruct the fermion operators, we must look for pairs of points of unit length 
in the weight lattice that differ by an element of the root lattice. If ai are the d simple 
roots of Spin(2d) and A* are the fundamental weights then: 

ai • = 6/, gij = ^di-dj, 5*-^' = 2A* • X^ (5.19) 

One also has the dual relation: 

{X%idi)g = 5pg (5.20) 

which will be important in what follows. Here p, q label the individual components of each 
vector, and the sum is over the vectors (not components). 

Eor the conformal dimensions, we have: 

= =^g"^ViVj (5.21) 
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Therefore to find operators of dimension {\ ■,\) we must look for sets of integers Vi for 
which: 

g^^ViVj = 2 (5.22) 

From Eq. ( |5.2C|) it follows that: 


g ^ {oti)p{otj')q — 2(5pq 


(5.23) 


hence the possible Vi are given by: 

v^f’ = {ai)p, p=l,2,---,(i (5.24) 


For the anti-holomorphic part we start by picking it independently from the same set: 

= {ai)q. However, we need to ensure that is \/2 times an integer. 

This is guaranteed by the fact that 

_ (5^.)^) = ^^((A^)^ - (A^)g) (5.25) 

which is indeed in the root lattice. We conclude that the fermion operators are given in 
bosonic language by: 

Dp{z)Dq{z) ^ Op^q = ( 5 . 26 ) 

where = \/2(A*)p. 

We can now look for the twist field, an operator cr^ satisfying Eq. 
a monodromy: 

(jfc : Dp{z) Dp{z) (5.27) 

corresponding to a shift: 

+ — (5.28) 

n 

(k) (k) 

This in turn will be induced by a shift (j)i ^ <pi + 27rQ ^ where Q ^ is a constant vector 
satisfying: 

u;(P)*cf) = ^ (5.29) 

for all p. Recalling that the last weight of Spin(2d) is A^*^) = (i ’ I ’''' > ^ easily hnd 
that the shift is given by: 

= ^(0>0>--- ,0,1) (5.30) 

Thus the twist field only acts on the last scalar Indeed, it takes the form: 


(2.2). This induces 




o 


Q{k)i 




(5.31) 
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where one must keep in mind that To check that this is correct, compute 

the OPE with Dp(z) to find: 




,(v)i/-(k), . k 

^ ^ = (z - Z )n 


(5.32) 


The crucial test of our twist field is whether it has the desired conformal dimension. We 
have: 


1 


{k)i^{k)j _ ijAk)Jk) _ ^ dd 


Using = i we get from which it follows that: 


fc=i 


E a ^ 

Ai. = — 
24 


k=-'- 


1 

n- 

n 


(5.33) 


(5.34) 


as desired. 

It remains to calculate the two-point function of each ak and thereby the replica 
partition function. Recall that the ordinary partition function for these theories is: 


Zi = 


1 


|^(T-)|2d 

1 1 




w^wGA-^y 

w—w^A^ 


E io> 


(5.35) 


12d 


2 I^WPEf2,3,4 


where, as usual, w'^ = gijw''w^ and similarly for w'^. The equivalence between the first 
and second lines of the above equation is the statement of Bose-Fermi equivalence between 
lattice bosons and multiple fermions with correlated spin structures. To demonstrate this, 
one starts with the observation that for Spin(2d), A^/ = Ar 0 Ay 0 A^ 0 Ac where 
^v,s,c are the lattices obtained by shifting by a fundamental weight in the vector, 
spinor and conjugate spinor representations respectively. Now we split a generic vector 
w G Aw into two classes: those that lie is Ar U Ay and those in As U Ac- In the former 
set we can write rc* = \/2 Xlp arbitrary integers Up, while in the latter set 

= [up + ^ )(A*)p where Up are again arbitrary integers. We similarly write re* = 

'^J2p'^pi^'')p and w'^ = V^Ylpi TUp 0 5 )(A*)p in the two respective cases. Finally, the 
restriction w — w £ A/j is implemented by inserting a projection operator that causes both 
w and w to lie in A^j or in Ay (in the first set) or both to lie in A^ or in Ac (in the second 
set). In this way we end up with the four Jacobi 0-functions 0,y(O|T), iz = 1,2, 3,4, of which 
01 vanishes for familiar reasons. 

For the un-normalised two-point function of twist fields, we find instead: 

. 2dfc2 


{{ak{z, z)a_k{0))) 


o'lm 

0i{i\r) 



w,w£Ay^ 

w—w^A^ 



(5.36) 
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Now we have 




(5.37) 


p=i 



It follows that: 




{{ak(z,z)a_k(0)}) 


1 e[{0\T) 

2 ^idlr) l??('r)p‘^ 


2d 


(5.38) 


Taking the product over k outside or inside the sum over spin structures gives us, respec¬ 
tively, the uncorrelated/correlated replica partition functions. As in the previous cases, 
under a modular transformation the replica partition function is covariant with the ex¬ 
pected prefactor. 

6. Conclusions 

The principle of modular invariance has been consistently useful in understanding confor¬ 
mal field theory. Here we have applied it to the study of Renyi entropies of free Dirac 
fermions and found that one can define more than one modular-covariant replica partition 
function. Our proposal is that the uncorrelated version, taking the product over repli¬ 
cas performed after the sum over spin structures, defines the correct R&yi/entanglement 
entropy as I" —?• 0 while the correlated version, taking the product over replicas before sum¬ 
ming over spin structures, defines the correct entropies as £ ^ L. We verified that in this 
way the holographically predicted relation to thermal entropy of Ref.|17|| is satisfied. Com¬ 
paring with results in the literature for free compact bosons, we showed that Bose-Fermi 
duality holds in the small and large-interval regimes. 

It is clearly of importance to obtain the correct modular-covariant replica partition 
function at arbitrary values of i. For a generic positive integer n one can write a variety of 
partition functions as follows. Consider the set S = {ki, k 2 , -'' j ^n} where ki = —1. 
Suppose that for a fixed value ki and a fixed spin structure u, we have a replica partition 
function Zn\ki). Partition S into subsets 5i, S 2 , • • • having mj elements of S in the jth 
set. Now we can define the replica partition function: 



( 6 . 1 ) 
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Clearly the fermion spin structures are correlated within each subset Sj but uncorrelated 
across different subsets. It is straightforward to verify that each of these is modular covari¬ 
ant, with a prefactor that can be eliminated as in Eq. ( ^.35 ). A particular case is when each 
subset Sj contains a single element kj. In this case we find the uncorrelated replica parti¬ 
tion function which we have argued is the correct one as £ ^ 0. Another particular case 
arises when there is a single subset, the entire set S. In this case one finds the correlated 
replica partition function which is correct as £ —>■ L. All other cases he somewhere in 
between. Therefore we conjecture that the full replica partition function of the free Dirac 
fermion is a linear combination (with ^-dependent coefficients) over all partitions: 


Zn — 


all partitions 


( 6 . 2 ) 


The physical intuition is that as the branch cut of length i expands from 0 to L, the 
degree of correlation among replicas steadily increases. The coefficients a must be such 
that only the terms corresponding to Z^ contribute at £ = 0, L respectively. It will be 
interesting to pursue this direction and try to determine the coefficients. 

This discussion strongly suggests that the modular-invariant replica partition function 
for free fermions is a very complicated quantity (in striking contrast to the one for a single 
spin-structure, which is relatively simple and very explicit for all i). From this point of 
view it is encouraging that the free boson result of Ref. [pl|] is also very complicated! After 
all, the two are supposed to be equal on putting ii = 1 in the latter. 

A possible application of the above ideas is to understand entanglement in more general 
conformal field theories, which at finite temperature and finite size remains largely an 
open problem. All the techniques used so far have been based on free field theory. In 
this context, it may be mentioned that minimal models and other CFT’s do have free- 
field representations with a screening charge |32| and their modular-invariant partition 
functions and correlators are quite well-known. It may be possible to determine replica 
partition functions for these theories, or at least their limiting behaviours for small and large 
intervals, using a combination of the free field representation and several physical principles 
including the relation to thermal entropy and the requirement of modular invariance. 
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Notation and conventions 

Our 2d theories are defined on a circular space of length L. The size of the interval over 
which we compute the Renyi entropy is denoted 1. The numbers Ri describe the radii of 
the scalar field variables. We use conventions where T-duality is implemented by R — 
and hence the self-dual radius is i? = \/2- The number of replicas is denoted n and used 
to calculate the nth Renyi entropy. The modular parameter is q = and the finite 
temperature case is described by r = where [3 is the inverse temperature. 
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